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Abstract 

To more fully take advantage of a low-cost, small footprint hybrid interfer¬ 
ometric/dispersive spectrometer, a mathematical reconstruction technique 
was developed to accurately capture the high-resolution and relative peak in¬ 
tensities from complex spectral patterns. A Fabry-Perot etalon was coupled 
to a Czerny-Turner spectrometer, leading to increased spectral resolution by 
more than an order of magnitude without the commensurate increase in spec¬ 
trometer size. Measurement of the industry standard Hg 313.1555/313.1844 
nm doublet yielded a ratio of 0.682, which agreed well with an independent 
measurement and literature values. The doublet separation (29 pm) is similar 
to the U isotope shift (25 pm) at 424.437 nm that is of interest to monitoring 
nuclear nonproliferation activities. Additionally, the technique was applied 
to LIBS measurement of the mineral cinnabar (HgS) and resulted in a ra¬ 
tio of 0.682. This reconstruction method could enable signihcantly smaller, 
portable high-resolution instruments with isotopic specihcity, benehting a 
variety of spectroscopic applications. 
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1. Introduction 


Numerous scientific, medical, industrial, and verification applications de¬ 
mand high-resolution spectral measurements, which can be accomplished 
using Fabry-Perot (FP) interferometry, among other methods [HEIE]- Re¬ 
cently, FP interferometry has been adopted in pulsed measurements via laser- 
induced breakdown spectroscopy (LIBS) [1]. Spectra from laser-produced 
plasmas are commonly measured by use of angular dispersion, such as in 
Czerny-Turner (CT) or echelle spectrometers [3]. In LIBS, high-resolution 
measurements are especially valuable for measurements of isotope ratios, 
where the typical isotopic shifts are on order a few to several tens of picome- 
ters. Currently, rapid and economical techniques with isotopic specificity 
are highly sought for nuclear safeguards and non-proliferation applications. 
Development of LIBS to detect and monitor facilities for illicit activities 
{e.g., isotope enrichment to undeclared levels) [HI E] is critical to risk pre¬ 
paredness and treaty verihcation. However, high-resolution isotope measure¬ 
ments presently require very large (few-meter scale) CT spectrographs p] 
or echelle spectrographs with complex instrument functions and limited dy¬ 
namic range [9]. More portable, simple, and cost-effective options are desir¬ 
able in held applications. A practical alternative is an inexpensive, compact 
FP etalon coupled to a much lower resolution CT imaging spectrometer |3]. 

Accurate analysis of the 2-D intensity pattern measured by the instrument 
is important for interpreting prohles of spectral lines. A variety of algorithms 
have been developed to deconvolve the input spectrum from the interferomet¬ 
ric measurement, but they rely on assumptions made about the shape {e.g., 
Gaussian, Lorentzian, Voigt functions) of the source spectrum (TUI EU 112] • 
Other techniques include analysis in the Fourier domain, which also assumes 
a known source distribution na. Direct least squares htting of 2-D ring im¬ 
ages from a single spectral line has been proposed, but is mainly aimed at 
calibrating and characterizing visible laser diodes [H] . Another alternative is 
data inversion [T3l[T3l[T7]. which can be effective in cases where the spectrum 
is unknown. Abbiss and Heeg measured spectra by varying the input angle of 
a tunable Liter [I6] . While such scanning methods are not suitable for pulsed, 
imaging LIBS applications, similar mathematical principles were used in this 
case. The reconstruction method employed to interpret the high-resolution 
spectrum captured by a hybrid interferometric/dispersive spectrometer takes 
into account the intensity distribution of fringes observed in a 2-D image, and 
can be applied to spectra of a broadband source without presuming the shape 
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of the input spectrum. In order to combine the benehts of 2-D algorithms 
with the speed necessary for rapid and high-resolution fringe measurement, 
this paper presents a simplihed and robust htting algorithm. The technique 
was also experimentally verihed with an actual FP fringe imaging system. 

The hybrid instrument, when used in conjunction with the reconstruc¬ 
tion technique discussed here, represents an effective tool for measuring the 
high-resolution spectrum with accurate peak intensity ratios. Rapid, high- 
resolution spectral reconstruction capability from this type of hybrid mea¬ 
surement is ideal in certain applications of LIBS. The resolution of the hy¬ 
brid interferometric-dispersive approach is at least an order of magnitude 
greater than the resolution of the CT spectrometer alone. This instrument 
can achieve the desired high resolution for measurement of the atomic iso¬ 
topic shift of uranium, while reducing the size and cost of the instrument by 
approximately an order of magnitude. 

2. Instrument model and spectral analysis 

The fringe pattern from an etalon exhibits angular and spectral depen¬ 
dence arising from constructive and destructive interference. A broad and 
complex spectrum incident on the FP etalon may produce a transmission 
pattern that cannot be uniquely reconstructed. This is due to degeneracy 
inherent in the interference pattern produced by the etalon, which exhibits a 
characteristic periodicity and wavelength-angle coupling, whereby an identi¬ 
cal phase shift can be produced by different combinations of wavelength and 
transmission angle. This degeneracy can be broken by coupling the etalon to 
a dispersive spectrometer, which provides an additional angular dispersion, 
as described by Effenberger et a/. |1]. A diagram illustrating the approach 
to image sampling and dispersion is shown in Fig. [^and the image incident 
on the 2-D detector array is shown in Figj^a). The high-resolution spectral 
information is encoded on the axis parallel to the spectrometer slit (he., or¬ 
thogonal to the dispersion plane of the spectrometer). The production of the 
interference ring pattern by the FP etalon can be described as 

poo 

B{e)= / T{e,x)A{e,x)dx, (i) 

Jo 

where A{6, A) is the spectrally and angularly dependent intensity of the light 
incident on the etalon and B{6) is the spectrally integrated intensity of the 


3 


transmitted light. The FP etalon transmission function is described by the 
Airy equation, 

T(6*, A) = (l + Fsin^ (27r(in/A cos 0)) (2) 

where F is the coefficient of hnesse, given by F = AR/ (1 — where R 
is the etalon reflectivity, d is the separation of the etalon mirrors, n is the 
refractive index in the space between the mirrors, A is the wavelength, and 
6 is the angle measured with respect to the etalon axis. In a typical FP 
optical arrangement, the transmission pattern is measured in the focal plane 
of a lens (far held). The fringe pattern measured in this way is comprised 
of concentric rings of light corresponding to the angularly and spectrally 
dependent interference maxima arising from propagation through the FP 
etalon. The maxima appear when the condition cos (6^) /A = /ctt/ { 27 rdn) is 
met (where k is an integer). Degeneracy arising from the coupling of angle 
and wavelength can be broken by introduction of angular dispersion, which 
leads to the measurement of 

B{e,x) = T{e,x)A{e,x). (3) 

Here, B{6, A) represents the recorded 2-D interference pattern. The measure¬ 
ment of this pattern is typically limited by the available spectral resolution, 
which is dictated by the properties of the arrangement used to produce dis¬ 
persion and the size of the detector pixels. At any wavelength Aq, the quantity 
measured is 

/^AoH-e/2 

mA)U=Ao=/ T{9,X)A{e,X)dX, (4) 

AAo-e/2 

where e is the spectral resolution at wavelength Aq. It is assumed for the 
purpose of this analysis that the resolution in the angular direction [6) is 
sufficient to accurately measure the shape of 5(6', A)|a=Ao- If is further as¬ 
sumed that light incident on the etalon does not exhibit angular dependence 
of the spectrum, so that the angular dependence of B{6, A) arises solely from 
the nonuniform angular intensity prohle and B{6) = T{6, A)A(A). The light 
is incident on the etalon from the output of an optical hber, and the Gaussian 
model chosen for the angular dependence of the intensity distribution is 

B{9) = Bo exp {-{9 - 9of/w^) + C, (5) 

where Bq is the peak intensity of the Gaussian envelope, 9o is the angle 
corresponding to the etalon axis, G is a constant, and w is the width of the 
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distribution envelope [IHl[in]. The envelope function can be used to correct 
the measured FP pattern for its angular intensity dependence by htting it to 
the Gaussian envelope B{9). 

This model of the measurement can be used for reconstructing the high- 
resolution features in the input spectrum directly from the measured pattern 
produced from the hybrid FP-dispersive instrument, without making addi¬ 
tional assumptions regarding the shape of the spectrum. In a real measure¬ 
ment, the discretization and normalization (to the best-fit envelope function 
A{6)) of the measured transmission pattern in the neighborhood of the wave¬ 
length Ao yields the vector 

Be = (6) 

Here, is the transmission matrix constructed by discretization of Eq. (|^, 
and Aa is the vector representing the spectrum in the vicinity of the wave¬ 
length Aq. The high-resolution spectrum can then be reconstructed as 

Aa = (7) 

which is, in practice, reduced to Ending the least-squares solution to this 
inverse problem. 

3. Experimental results and discussion 

The reconstruction method was applied to data collected from a hy¬ 
brid FP instrument that used a 550 mm focal length imaging spectrometer 
(Horiba Jobin Yvon) equipped with a 1800 mm“^ grating and a 1024x1024 
pixel array iStar ICCD (Andor) composed of 13 pm pixels. Coupled to 
the entrance of the spectrometer was a FP etalon (SLS Optics). Light was 
transported to the etalon by a multimode 400 pm core diameter optical fiber 
(Ocean Optics). Fiber-optic connector, etalon, and lens were affixed in a cage 
setup, with multi-axis adjustments to enable accurate alignment of the cir¬ 
cular transmission pattern to the spectrometer slit. A 25 mm diameter lens 
with a focal length of 100 mm was used to focus the output of the etalon onto 
the entrance slit of the spectrometer. The light source used was a Hg(Ar) 
calibration pen lamp (Oriel). 

The FP etalon was coupled to the entrance slit of the spectrometer to pro¬ 
duce fringes, as shown in Fig. [^a). The section of the spectrum correspond¬ 
ing to the 313.1555/313.1844 nm doublet was summed along the horizontal 
axis to obtain the 1-D profile of the transmitted intensity, shown in Fig. |^b). 
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The symmetric pattern was collapsed to summing about the etalon axis to 
improve the signal-to-noise ratio (SNR). A Gaussian envelope function was 
fit to the fringes exhibiting a characteristic double-peaked shape to obtain the 
A{9) correction, as shown in Fig. [^c). The normalized transmitted intensity 
pattern is shown in Fig. |^d). 

It is necessary to calibrate the angular axis before formulating the in¬ 
strument function. This calibration was performed using the neighboring 
312.5674 nm line (Fig.j^a)). The angular position of highest intensity (peak) 
in the fringe profile was used to perform the calibration in Fig. Because 
the paraxial approximation is valid, the calibration can be described by a 
linear function, also shown in Fig. The region of the fringes selected for 
this £t did not include the area near the optical axis or the outermost fringes, 
which may exhibit distortions due to their low intensity and possible optical 
aberrations. 

Next, the instrument function for the 313 nm doublet was determined. 
The ranges of 6 and A used to construct the instrument function were 20 mrad 
< 9 < 32 mrad and 313.106 nm < A < 313.234 nm, respectively. The etalon 
had air-spaced (n = 1) mirrors separated by a distance oi d = 0.44 mm. 
The etalon reflectivity {R = 0.73) was determined by fitting the measured 
intensity prohle to the known FP etalon transmission function, Eq. ([^, as 
shown in Fig. 

Interpolation between points and increased sampling of the experimental 
spectrum was performed prior to reconstruction of the source spectrum, Aa. 
The least-squares solution to the system of linear equations associated with 
Eq. ([^ could have been used to reconstruct Aa. However, Eq. ([^ is an 
example of a problem where a is an ill-conditioned matrix HT]. Hence, 
the pseudoinverse of the matrix was used to calculate a least-squares 
solution to the system of equations using 

= Aa. (8) 

The pseudoinverse solution was implemented in Mathematica (Wolfram Re¬ 
search) using the function pseudoinverse. Because of the ill-condition of 
the matrix T^^a, it was necessary to specify the tolerance parameter. The 
singular values matrix are the non-negative square roots of the eigenvalues 
of Tg aT^.a, where Tg is the conjugated transposed matrix. Singular values 
smaller than the product of the tolerance and the largest singular value of 

A were omitted [20]. In practice, the choice of tolerance controls the influ¬ 
ence of smaller singular values in the composite error of the reconstruction. 
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and is a trade-off between redncing ronnd-off error at the expense of a larger 
residnal. A tolerance of 0.1 was nsed thronghont the analysis to yield the 
optimal approximation to the solntion. 

Minimization of the residnal error of the reconstrnction was the next 
step. The widths of Hg lines from a low-pressnre discharge calibration lamp 
are approximately 1 pm [5]. The width of the reconstrncted spectrnm is 
dictated primarily by the reflectivity R. Becanse the widths are narrow 
compared to the instrnment fnnction, an accnrate reconstruction assuming 
this experimentally determined reflectivity was not possible. By assuming 
a broader source spectrum, it was possible to minimize residual error in 
the spectral reconstruction. A reflectivity value of = ai?, where a is 
a multiplicative constant of R that minimizes the residual error, || B^i — 
Ta,6»Aa II was assumed. An optimal value a = 1.16, obtained by scanning a 
range of values 1.00< a <1.30, was utilized to minimize the residual error of 
the reconstruction. 

A spectrum of Hg measured using the CT spectrometer without the FP 
etalon is shown in Fig. 4(a) for comparison. In this case the resolution 
(~70 pm) was not sufficient for resolving the 313 nm doublet. The high- 
resolution reconstruction provided by the hybrid instrument is shown in 
Fig. Qb) and summarized in Table 1. The resolution of the 313.1844 nm 
line was ~4 pm, which exceeds the resolution of the CT spectrometer by 
more than one order of magnitude. These results demonstrate that, while 
the resolution results from the design of the entire hybrid instrument, the res¬ 
olution enhancement is a consequence of the introduction of the FP etalon. 
An accurate reconstruction method is critical to fully realize the utility of a 
hybrid instrument. For comparison, the 313.1844 nm line in Fig. |^d) has 
a resolution of ~8 pm using more traditional methods for converting the 
fringe pattern to a spectrum. Additionally, Effenberger et al. estimated the 
resolution to be ~11 pm from the measured peak separation in spectra ac¬ 
quired with an Hg(Ar) lamp coupled into a FP attached to a 500 mm CT 
spectrometer with a 512x512, 24 pm/pixel detector array. 

Correcting for the angular drop in intensity of the measured FP pat¬ 
tern was critical for peak ratio determination; the reconstructed peaks of 
non-normalized data yielded a ratio of 0.629, which is lower due to the de¬ 
creasing intensity as a function of distance from the interferogram center. 
Validation of the reconstructed hybrid instrument spectrum was performed 
by comparing the results with those acquired with a high-resolution echelle 
spectrometer (DEMON, Lasertechnik Berlin) using the same Hg(Ar) lamp. 
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The DEMON results are also provided in Fig. and Table along with 
literature values m 

Averages reported in Table [T] are based on 10, 5, and 20 spectra for 
sources 1, 2, and 3, respectively. The reconstructed peak wavelengths were 
313.156 nm and 313.185 nm for the Hg doublet, and agree well with actual 
values. While the resolution is not as high for the hybrid instrument, the 
spectra are similar (Fig.^b)) and give a very similar ratio accuracy (Table [^. 


Table 1: Peak ratios of the 313 nm doublet of a Hg(Ar) calibration lamp measured by 
different sources. The same lamp was measured by method 1 and 2. 


Source 

Ratio 

Uncert.(%) 

A/AA 

1. 

Hybrid device 

0.682 

1.8 

O 

O 

o 

CO 

2. 

DEMON 

0.681 

2.4 

-75,000 

3. 

Ref. [21] 

0.687 

7.2 

N/A 


To further explore the utility of this spectral reconstruction method, it 
was applied to the analysis of LIBS data from the mineral cinnabar (HgS) 
acquired with a different hybrid spectrometer |3]. In this case, plasma was 
generated by 1064 nm laser pulses at 25 mJ per pulse. The spectrum was 
collected from the accumulation of 600 laser pulses at a rate of 10 Hz and 
under 10 Torr He atmosphere, which is the optimum pressure for the Hg 
313 nm doublet [1]. The reconstructed doublet, shown in Fig. |^b), exhibits 
a ratio of 0.682 and is in good agreement with values reported in Table 
Peak broadening in the LIBS spectrum is due to plasma dynamics |22j . 

4. Conclusion 

This spectral reconstruction approach for hybrid FP/dispersive spectrom¬ 
eters demonstrates that simple, low-cost instruments with high resolution and 
accurate relative peak intensities can be realized. Accurate wavelength and 
ratio determination of Hg lines using this instrument has been presented, 
where the accuracy depends on proper normalization and calibration of the 
measured interferogram. Because only a fraction of light is coupled to the 
CT spectrometer, a fraction of the signal is lost in this method. Summation 
of signal from different orders produced by the FP etalon may be able to 







increase the SNR, but was not implemented in the reconstruction discussed 
here. Experimentally, a higher repetition rate laser may be used to speed 
collection time and further increase the SNR. While such instruments would 
be invaluable for LIBS measurements requiring accurate isotope ratio mea¬ 
surements, the technique could enable practical instruments for a variety of 
high resolution held applications. 


9 



The authors would like to thank Andrew J. Effenberger, Jr. and Elizabeth 
J. Judge for their assistance. Research was performed under appointment to 
the Nuclear Nonproliferation International Safeguards Graduate Fellowship 
Program sponsored by the National Nuclear Security Administration’s Next 
Generation Safeguards Initiative (NGSI). Material is based upon work sup¬ 
ported by the U.S. Department of Homeland Security under Grant Award 
Number, 2012-DN-130-NF0001-02 and in part funded by the Gonsortium for 
Verihcation Technology under Department of Energy National Nuclear Secu¬ 
rity Administration award number de-na0002534. The views and conclusions 
contained in this document are those of the authors and should not be in¬ 
terpreted as necessarily representing the official policies, either expressed or 
implied, of the U.S. Department of Homeland Security. Research was also 
sponsored by the U.S. Department of Energy under DOE Idaho Operations 
Office Gontract DE-AG07-05ID14517. 


10 



References 


[1] P. B. van der Wal, J. Slingerland, High resolution Fabry-Perot spec¬ 
trometer for infrared astronomy, Appl. Opt. 18 (1979) 641-648. 

[2] M. D. Zhou, C. Yang, Z. Liu, J. P. Cysyk, S. Y. Zheng, An implantable 
Fabry-Perot pressure sensor fabricated on left ventricular assist device 
for heart failure, Biomed Microdevices 14 (2012) 235-245. 

[3] D. R. Vogus, V. Mansard, M. V. Rappa, T. M. Squires, Measuring 
concentration helds in microfluidic channels in situ with a FabryPerot 
interferometer. Lab Chip 15 (2015) 1689-1696. 

[4] A. J. Effenberger Jr., J. R. Scott, Practical high-resolution detection 
method for laser-induced breakdown spectroscopy, Appl. Opt. 51 (2012) 
B165-B170. 

[5] C. Aragon, J. A. Aquilera, Characterization of laser induced plasmas by 
optical emission spectroscopy: A review of experiments and methods, 
Spectrochim Acta B. 63 (2008) 893. 

[6] C. Hanson, S. Phongikaroon. J. R. Scott, Temperature effect on laser- 
induced breakdown spectroscopy spectra of molten and solid salts, Spec¬ 
trochim. Acta B. 97 (2014) 79-85. 

[7] P. Ko, K. C. Hartig, J. P. McNutt, R. B. D. Schur, T. W. Jacomb- 
Hood, 1. Jovanovic, Adaptive femtosecond laser-induced breakdown 
spectroscopy of uranium. Rev. Sci. Instrum. 84 (2013) 013104. 

[8] C. A. Smith, M. A. Martinez, D. K. Veirs, D. A. Cremers, Pu-239/Pu- 
240 isotope ratios determined using high resolution emission spec¬ 
troscopy in a laser-induced plasma, Spectrochim. Acta B. 57 (2002) 
929-937. 

[9] D. A. Cremers, A. Beddingheld, R. Smithwick, R. C. Chinni, C. R 
Jones, B. Beardsley, L. Karch, Monitoring Uranium, Hydrogen, and 
Lithium and Their Isotopes Using a Compact Laser-Induced Breakdown 
Spectroscopy (LIBS) Probe and High-Resolution Spectrometer, Appl. 
Spectrosc. 66 (2012) 250-261. 


11 



[10] H. P. Larson, K. L. Andrew, A least squares deconvolution technique 
for the photoelectric Fabry-Perot Spectrometer, Appl. Opt. 6 (1967) 
1701-1705. 

[11] J. Koo, J. Kim, Deconvolution algorithm for a Fabry-Perot interferom¬ 
eter, App. Opt. 26 (1987) 1094-1099. 

[12] D. A. Naylor, N. Lermer, I. Furniss, Deconvolution of Fabry-Perot spec¬ 
tra, Infrared. Phys. 31 (1991) 401-408. 

[13] V. G. Cooper, Analysis of Fabry-Perot Interferograms by means of their 
Fourier transforms, App. Opt. 10 (1971) 525-530. 

[14] M. OHora, B. Bowe, and V. Toal, Least squares algorithm for rapid 
analysis of Fabry-Perot fringe patterns, J. Opt. A: Pure Appl. Opt. 7 
(2005) S364-S367. 

[15] D. M. Zipoy, Fabry-Perot inversion algorithm, App. Opt. 18 (1979) 1988. 

[16] B. Heeg, and J. B. Abbiss, Piezospectroscopic imaging with a tunable 
Fabry-Perot Liter and Tikhonov reconstruction. Opt. Lett. 7 (2007) 
859-861. 

[17] J. B. Abbiss, B. Heeg, Imaging piezospectroscopy. Rev. Sci. Instrum. 79 
(2008) 123105. 

[18] T. J. Scholl. S. J. Rehse, R. A. Holt, S. D. Rosner, Broadband precision 
wavelength meter based on a stepping FabryPeerot interferometer. Rev. 
Sci. Instrum. 75 (2004) 3318. 

[19] M. C. Hirschberger, G. Ehret, Simulation and high-precision wavelength 
determination of noisy 2D Fabry-Perot interferometric rings for direct- 
detection Doppler lidar and laser spectroscopy, App. Phys. B 103 (2011) 
207-222. 

[20] G. H. Golub, C. Reinsch, Singular value decomposition and least squares 
solutions, Numer. Math. 14 (1970) 403-420. 

[21] J. Reader, C. J. Sansonetti, J.M. Bridges, Irradiances of spectral lines 
in mercury pencil lamps, Appl. Opt. 35 (1996) 78-83. 


12 



[22] J. R. Scott, A. J. Effenberger Jr., J. J. Hatch, Laser-induced Breakdown 
Spectroscopy-Theory and Applications (Springer-Verlag, 2014), Chap. 

4. 


13 



Optical Etalon 



spectrometer 


Figure 1: Schematic diagram of optical detection setup. 
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Figure 2: (a) Measured pattern corresponding to the (I) 312.5674 nm line and (II) 
313.1555/(111) 313.1844 nm doublet; (b) 1-D integrated profile of doublet; (c) intensity en¬ 
velope function for (II) 313.155 nm (dots - experimental data; dashed line - ht to Eq. (§); 
(d) normalized 1-D proHle. 


15 


































Figure 3: Best-fit (solid line) of the measured 312.5674 nm first-order interference peak by 
Eq. (§ yields the reflectivity of the etalon {R = 0.73). Experimentally measured points 
are shown as black markers. Inset shows calibration of the angular (9) axis using the 
location of four 312.5674 nm constructive interference peaks (black markers), The linear 
fit is shown as a dashed line. 
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Figure 4: (Color online) (a) The (I) 312.5674 nm line and (II) 313.1555/(III) 313.1844 nm 
doublet measured with the CT spectrometer alone; the doublet (II+III) was not resolved, 
(b) The spectrum of the (II) 313.1555/(111) 313.1844 nm doublet (black, dotted) recon¬ 
structed from the hybrid FP/CT measurement compared to DEMON measurement (blue, 
solid). LIBS of HgS reconstructed from a hybrid FP/CT measurement is also shown (red, 
dashed). 
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